EPJ manuscript No. 

(will be inserted by the editor) 



Bimodalities: a survey of experimental data and models 



o 
o 

> 
O 

:^ 

o 



O. Lopez^ and M. F. Rivet^ 

^ Laboratoire de Physique Corpusculaire, IN2P3-CNRS/ENSICAEN /Universite, F- 14050 Caen cedex, France 
^ Institut de Physique Nucleaire, IN2P3-CNRS, F-91406 Orsay cedex, France 

February 8, 2008 

Abstract. Bimodal distributions of some chosen variables measured in nuclear collisions were recently 
proposed as a non ambiguous signature of a first order phase transition in nuclei. This section presents a 
compilation of both theoretical and experimental studies on bimodalities performed so far, in relation with 
the liquid-gas phase transition in nuclear matter. 
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After a formulation of the theoretical bases of bimodal- 
ity, world-wide experimental results will be reviewed and 
discussed, as well as the occurrence of some kind of bi- 
modality in models. Finally conclusions on the perspec- 
tives of such analyses in the near future and the possi- 
ble connections to other proposed signals of the liquid-gas 
phase transition in nuclear matter will be given. 



1 Theoretical bases 
1.1 Definition 

Bimodality is a property of finite systems undergoing a 
first order phase transition [Tl[2l[3] . It is thus a generic fea- 
ture which concerns not only Nuclear Physics but a broad 
domain of physics such as Astrophysics, or Soft Matter 
Physics. Bimodality means that the probability distribu- 
tion of an order parameter of the considered system at 
phase transition exhibits two peaks separated by a mini- 
mum. Indeed, if the system is in a pure phase, the order 
parameter distribution consists in one peak and can be 
characterized by its mean value and its variance. By con- 
trast, if the system is in the coexistence region, the dis- 
tribution presents two peaks, well separated, whose prop- 
erties are related to the two different phases of the sys- 
tem [5]. Bimodality is then one of the signals associated 
to a first order phase transition [3] , beside others such as 
scaling laws, critical exponents or negative heat capacities. 

In the following, the term "bimodality" will abbreviate 
"the probability distribution of some variable, in a given 
region of the phase diagram of the system, is bimodal" . 



1.2 Pioneering studies 

Bimodality and its relationship to phase transition has 
been studied since the 80's. Fig. [1] shows an Ising model 
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Fig. 1. Evolution of the magnetization M, as a function of 
the applied magnetic field H , in the Ising Model for a lattice 
defined by the size L. The bottom panel presents a schematic 
probability distribution of the magnetization between —Ml 
and -\-Ml around the critical field value He- Taken from [4]. 
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simulation of a fcrromagnct studied by Binder and Lan- 
dau [4]. In this analysis, the authors studied the mag- 
netization M of the system as a function of the apphed 
magnetic field H. When the magnetic field comes close 
to the critical value ffc, the spontaneous magnetization 
of the ferromagnet presents a sudden change; in this case 
the probability distribution of the magnetization is never 
bimodal, as the system "jumps" suddenly from the nega- 
tive value -Mgp to the positive one +Msp'. the transition 
between the two regimes is sharp at the thcrmodynamical 
limit. 

By contrast, when the size of the system is finite (and 
defined by the number of sites L), the step function is 
replaced by a smooth curve in fig[l] with a slope propor- 
tional to - where d is the dimensionality of the system. 
Consequently, in the vicinity of He, the magnetization M 
exhibits a bimodal structure as shown in the bottom panel 

of fig. m 

1.3 Link with phase transition in Thermodynamics 

It was recently demonstrated by Chomaz and GulmincUi 
that bimodality of the probability distribution of the order 
parameter is equivalent to the other definitions of phase 
transition proposed up to now [5]. 

Relationship to the Yang-Lee theorem The Yang-Lee 
theorem [6] is considered as the standard definition of 
first-order phase transitions at the thermodynamic limit. 
As demonstrated in [5] bimodality is a necessary and suf- 
ficient condition for zeroes of the partition sum in the 
control intensive variable complex plane to be distributed 
on a line perpendicular to the real axis. 

Anomaly of thermodynamical potentials A first order 
phase transition is characterized by an inverted curvature 
of the relevant thermodynamical potential (entropy, free 
energy) [TJIH] . This feature is also equivalent to a bimodal- 
ity in the event probability of the given order parameter 
X as displayed in the left part of fig. [H 

Negative derivatives of the thermodynamical potentials 

A first order phase transition was also related to a back- 
bending in the equation of state of the system [7] , charac- 
terized by a negative second derivative of the thermody- 
namical potential, as for example the heat capacity if the 
energy is the order parameter, fig. [2] right. 

1.4 Microcanonical vs. canonical ensemble 

Among the observables signing a phase transition, the 
heat capacity is related to the fluctuations of the par- 
tial energy of the system and need to be studied in the 
microcanonical ensemble, while bimodality can only be 
observed when the system is free to fluctuate in terms of 
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Fig. 3. Energy distributions obtained in the Gaussian ensem- 
ble for different a = N/N'. Taken from |9]. 

the associated extensive variable (i.e. energy or volume). 
This case corresponds to canonical or isobar ensembles. In 
other words, events must be selected without constraint 
on the extensive variable in order to study bimodality. 
However in nuclear physics experiments, the two colliding 
nuclei form an isolated system: it seems thus natural to 
work in a microcanonical ensemble, and cuts can be ap- 
plied on the energy of the system, determined for instance 
by calorimetry. It seems conversely out of reach to be in 
a canonical framework which would require the existence 
of a large heat bath. 

The situation is not hopeless, as it was shown some 
years ago that properties of phase transitions can be ob- 
served even if the working ensemble is not strictly micro- 
canonical or canonical, but is an interpolating ensemble. 
In Gaussian ensembles for instance, it is supposed that 
N particles are in contact with a system of N' particles 
acting as a heat bath at temperature T. When N' varies 
from to oo, the working ensemble mimicks the transition 
between microcanonical and canonical [5]. fig. [3] presents 
the results of such a simulation, where it is clearly seen 
that the probability distribution of the energy - in the 
transition region - presents a bimodal shape only when 
N/N' is small enough (< 1/1000) while for larger N/N', 
the situation is that of the microcanonical case with only 
one peak in the distribution. 



1.5 Liquid-Gas phase transition 

Since nuclei are supposed to undergo a liquid-gas phase 
transition, specific studies of this peculiar transition were 
undertaken through Lattice-Gas calculations. In liquid- 
gas phase transitions, volume as well as energy are order 
parameters. The bimodality of the event probability dis- 
tribution in the first order phase transition region is ev- 
ident in fig. m which shows the location of events in the 
Volume vs. Energy plane (top left). The projections along 
the axes {E, V) also display the expected bimodality, as 
does a linear combination of these two order parameters 
(bottom right, red curve). In this framework (Lattice-Gas 
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Fig. 4. Probability distributions of the energy E, volume V 
and a combination of the two variables coming from a Lattice- 
Gas simulation in the canonical ensemble. Taken from [2]. 



model), bimodality is evidenced if we are able to select 
(sort) events in a canonical way (or as close as possible, 
see previous section), and plot the event probability distri- 
bution of the energy or volume, or any observable directly 
related to them. 



2 Experimental observations 

Since bimodality was proposed as a signature of liquid- 
gas phase transition, it was extensively searched for in 
event samples resulting from nuclear collisions; studies 
were made for central collisions, where the liquid-gas phase 
transition is clearly evidenced by previous analyses (see 
chapter "phase transition") as well as for peripheral col- 
lisions, where a large range of excitation energy can be 
explored. 



2.1 Central collisions: systems with mass ^ 250 

Systems with total mass close to 250 were studied with the 
INDRA array using two entrance channels, an asymmetric 
one, Ni-f Au, and an almost symmetric one, Xe-|-Sn. In 
both cases, in the incident energy range scanned, it was 
shown that a fused system was formed in central collisions. 

Bellaize et al |lU| have reported the observation of bi- 
modality of the size asymmetry of the two largest frag- 
ments in central events for the Ni-|-Au system at 32A, 52A 
and 90AMeV. It was associated with two fragmentation 
patterns (see first row of fig. [5]), one similar to residue- 
evaporation (one large fragment with few small ones, zone 1 
in fig. [5]), the other to multifragmentation (fragments of 
nearly equal size, zone 2). A variable built with the charges 
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Fig. 6. Probability distributions of the charge asymmetry be- 
tween light (Z=3-12) and heavy fragments (Z>12) for fused 
events in the Xe-|-Sn system at 32yl, 39A, 45A and SOvlMeV. 
Taken from [11] ■ 
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Fig. 7. Charge asymmetry obtained by using a Stochastic 
Mean-Field simulation (BOB [12]) for central events of the 
Xe-|-Sn system at 32A MeV. Unpublished results from the au- 
thors of [H. 



of the three largest fragments, Zi, Z2, Z3 in decreasing or- 
der, 

Zi - 3(^2 + Z3) (1) 

also has a bimodal distribution at 32A and 52AMeV, 
as shown in the bottom row of fig. (S] but no longer at 
90A MeV. This fact is compatible with the location of the 
system in the coexistence region below 52AMeV, where it 
can experience a first order phase transition by exploring 
different densities and temperatures. For higher energies 
(here 90 A MeV), the system passes directly through the 
coexistence region and we observe only the presence of 
the multifragmentation regime, which could indicate that 
the system explores only the low density part of the phase 
diagram. 

fig.|6]shows the distributions obtained when looking at 
the asymmetry ratio between heavy, [Z > 13), and light, 
(Z=3 — 12), fragments 



(2) 



for single source events produced in central Xe+Sn colli- 
sions between 32 A and 50^ MeV [ll] . Bimodality is present 
at all energies, with dominant "liquid-type" events at 
32AMeV, and a dominance of "gas-like" events at and 
above 45AMeV; the two types of events are in roughly 
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Fig. 5. Correlation between the two largest fragments, Z\ and ^2 obtained in central collisions for the Ni+Au system at J,2A 
(left), 52^4 (middle) and 90^MeV (right). The bottom row shows an asymmetry variable built as a linear combination of the 
atomic number of the three largest fragments. Taken from [10) . 



equal number at 39AMeV, where other phase transition 
signals have been already observed (see chapter "Many 
fragment correlations"). The authors of [TT] relate the cho- 
sen asymmetry variable to the density difference between 
the coexisting liquid and gas phases of nuclear matter. 
The same variable was built for the events resulting from 
a Stochastic Mean Field simulation [12] of head-on col- 
lisions between Xe and Sn at 32^MeV. In this simula- 
tion, which was shown to well reproduce many experimen- 
tal features, single variable distributions as well as differ- 
ent correlations [T51[nifT5] (see chapter "Many fragment 
correlations" ) , the system enters the coexistence region 
and multifragments through spinodal decomposition. The 
equivalent of fig. [S] for simulated events is shown in fig. [71 
the picture is very similar to the experimental data at the 
same energy (black stars in fig. [5]) , a bimodal behaviour 
appears with a dominance of events of liquid type. 



2.2 Central collisions: systems with mass ^ 100 

Central collisions between two ^^Ni nuclei were studied 
at incident energies between 32 and 90AMeV; event se- 
lection was made through a discriminant factorial analy- 
sis trained, at variance with rcf [16] . on the complete ex- 
perimental events. A bimodal distribution of the largest 
fragment was observed at 52^ MeV, intermediate between 
the Gaussian distributions measured at lower energies and 
the asymmetric distributions found from 74A MeV up [T7] . 
The minimum is rather shallow (about 80% of the peak 
value); at 64ylMcV a bimodal distribution persists, but 
now the peak on the more fragmented side is dominant. 
Conversely the distributions of the fragments of higher 
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Fig. 8. Distributions of the largest fragment for central Ni-|-Ni 
collisions from 32 to 90^4 MeV (bottom). The same distribu- 
tions at the four lowest energies are displayed in linear scale in 
the top panels. Taken from |17] . 



rank (not shown) are monotonous. To our knowledge, it 
is the only direct observation of bimodality on the largest 
fragment. 
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2.2.1 Going further 

Central collisions allowed to study and evidence a bimodal 
behaviour of some asymmetry variables, which can be con- 
nected to the density difference between a liquid and a 
gas phase; in that sense they would be good candidates 
for being order parameters of a liquid-gas type transition. 
Nevertheless, several drawbacks can be pointed out; firstly 
it was shown that the lighter fragments exhibit a pre- 
equilibrium component in Ni4-Au [10] , while radial flow 
effects were recognised in symmetric systems, Xe-|-Sn [T8l 
fT9ll20] and Ni-^Ni [17]. But above all, the sorting of cen- 
tral events selects a rather narrow region in excitation 
energy for each incident energy (about 1-2A MeV at half 
maximum of the distribution). This is closer to a micro- 
canonical working ensemble and may prevent a very clear 
observation of bimodality. 

2.3 Quasi-projectiles in peripheral collisions 

Analyses of quasi-projectilcs formed in peripheral and semi- 
peripheral reactions arc thus mandatory, as they allow to 
overcome some of the abovementioned problems. In par- 
ticular a broad excitation energy distribution of quasi- 
projectilcs (QP) can be accessed. Exchanges of energy and 
particles with the quasi-target (QT), while it lies in the 
neighbourhood of the QP and especially when it is heavy, 
mimick a small heat bath and an almost canonical sorting 
can be envisaged. Whenever the incident energy is high 
enough, the different components (the QT and the QP, 
and the pre-cquilibrium or neck part) can be better dis- 
entangled, or at least the uncertainties caused by their 
existence can be circumvented. 

Most of the studies on quasi-projectiles arise from Au 
on Au collisions at various energies. Extensive results con- 
cerning a very light nucleus, close to Argon were also 
recently proposed. Several variables are used for sorting 
events as a function of the violence of the collisions; among 
the most commonly employed one can cite multiplicities 
and the transverse energy (relative to the beam axis) of 
charged products, either all of them or only light charged 
particles (Z = 1, 2) [21,22,23,24. Other sorting are based 
on Zbound (the sum of charges for fragments, Z > 2), as 
proposed by the ALADIN collaboration [25j, or on the 
excitation energy (NIMROD collaboration) [26j . 

2.3.1 Au quasi-projectiles at relativistic energies 

The ALADIN collaboration reported the presence of bi- 
modality for peripheral Au-|-Au reactions at lA GeV [27] . 
Fig. [9] shows the Zbound distribution (top panel) where 
is highlighted the selected region, Zbound = 53 — 55, for 
which was drawn the charge asymmetry between the three 
largest fragments 

Z1-Z2- Z3 (3) 

in the bottom panel. The charge asymmetry exhibits two 
components, the first one centered at low values (close to 
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Fig. 9. Tihound (top) and charge asymmetry distributions (bot- 
tom) for the Au+Au system at lA GeV. The bottom panel cor- 
responds to the Tibound selection displayed by the highlighted 
area in the top panel. Taken from [27] . 

0), which is associated to multifragmentation events, and 
the second one located at values around 40, which is more 
likely due to an evaporation residue of charge Z close to 
Zbound- It is worth saying that a percolation simulation 
was able to reproduce this bimodality in the charge asym- 
metry at the transition point. In this case, this is a second 
order phase transition. This point will be discussed in the 
"pending questions" section. 



2.3.2 A smaller system with mass ^ 40 

In a very complete analysis. Ma et al |26| scrutinized data 
collected with the NIMROD array. They were able to re- 
construct, from their emitted particles and fragments, the 
quasi-projectilcs formed in 47^1 MeV Ar-|-Al,Ti,Ni colli- 
sions. The method used consisted in tagging the particles 
with the help of a three moving source fit (QP, QT and 
mid-rapidity) and then attributing to each of them, event 
per event, a probability to be emitted by one of these 
sources. Completeness of quasi-projectiles, {Zqp >12), 
from semi-peripheral collisions was further required; QP 
excitation energy was determined using the energy bal- 
ance equation. The distributions of excitation energy so 
obtained for the three targets superimpose, showing that 
the QP excitation energy calculation is under control. 

Plots of the charge of the second largest fragment vs. 
the largest one are shown in fig. (TUl As for heavier sys- 
tems, the topology evolves from residue-evaporation to 
multifragmentation with increasing excitation energy. An 
equipartition of events between two topologies is observed 
for E* /A =5.5 MeV, where at the same time fluctuations 
on the size of the largest fragment are the largest, the 
power-law exponent for the charge distribution is mini- 
mum, and scaling laws are present. Here again, bimodality 
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Fig. 10. Correlation between the largest charge (Zmax) and the second largest {Z2ndmax) for the QP in peripheral Ar+Ni 
collisions at 47AMeV. Panels from Excl to Exc9 correspond to a selection in increasing excitation energy (see text). Taken 
from [26]. 



is observed at the same time as other possible indicators 
of a phase transition. 

2.3.3 Toward a canonical event sorting? 

In the previous cases the sorting for peripheral reactions 
uses properties of the studied source itself (here the QP) 
and is then probably more akin to a microcanonical than a 
canonical sorting. Indeed the bimodal character of the dis- 
tribution is not very marked, as expected if the experimen- 
tal sorting constrains strongly the excitation energy p5] . 
To attempt a true canonical sorting one must discriminate 
the studied system from some heat bath. A first tentative 
in that aim was the study of Au quasi-projectiles through 
a sorting performed on the transverse energy of the par- 
ticles of the Au quasi-target (as the system is symmetric, 
this amounts to particles emitted backward in the cm.). 
This sorting is illustrated by the results presented here- 
after. 

2.3.4 Au-like nuclei in a "canonical" sorting 

Au quasi-projectiles from Au-|-Au collisions at various in- 
cident energies were widely studied. Two examples are 
given here, at 35AMcV - results from the MULTICS- 
MINIBALL collaboration [29] - and at SOAMeV, data 
from the INDRA/ALADIN collaboration [3D]. In both cases 
data were sorted versus the transverse energy of the QT 
light charged particles. The charge of the largest fragment 
in each event is plotted in figs I11I12I versus the charge 
asymmetry of the two largest fragments, 

(Zi-Z2)/(Zi+Z2). (4) 

Whatever the incident energy the picture evolves from 
an evaporation residue to a multifragmentation configu- 
ration, passing through a zone where the two topologies 



coexist, separated by a neat minimum; in this zone (last 
one at 35 A MeV, third one at 80^ MeV) the distributions 
present a bimodal behaviour. Note that the bimodal char- 
acter is not very strong when one projects the bidimen- 
sional figures on either Z^aa or on the asymmetry. This 
is attributed in |30| to the presence of preequilibrium ef- 
fects, and some remaining aligned momentum which tend 
to shallow the minimum of a bimodal distribution. 



3 Bimodality in models 

Different statistical as well as dynamical models explicitely 
or implicitely contain a phase transition. They predict the 
occurrence of bimodal distributions for selected variables 
around some transition energy. Examples are given in this 
section. 



3.1 SMM: Statistical Multifragmentation Model 

Buyukcizmcci, Ogul and Botvina analyzed SMM sim- 
ulations for heavy nuclei of various sizes, with excitation 
energy ranging from 2 to 20 MeV/nucleon. They found 
that all nuclei exhibit the same caloric curve, depicted in 
the top panel of fig. [131 with the well-known "plateau" 
between 4 and 7 MeV/nucleon (note in passing that the 
common temperature at plateau whatever the mass of the 
considered nucleus is in contradiction with the experimen- 
tal results analyzed in ref [32|). In the same energy interval 
as the plateau, the fluctuations of A^ax (not shown) and 
of the temperature (panel (b) of the figure) are maximum. 
The authors sorted the events following the size of the 
largest fragment, Amax- They defined two event classes, 
one with A^ax > 2Ao/3, representative of the residue- 
evaporation channel and the other one with A„iax < Ao/3, 
characterizing multifragmentation events - Aq being the 
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Fig. 11. Correlation between the charge of the largest frag- 
ment (Zi) and the charge asymmetry (asymi2) between the 
two largest fragments for peripheral events of the Au+Au 
system at 35^4 MeV. The panels correspond to a selection 
in increasing transverse energy of particles coming from the 
QT side from top left to bottom right. Taken from [29j . 
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Fig. 12. Same as for fig. [TT] but for Au+Au at SOAMeV. 
The panels corresponds to the same cuts in QT transverse 
energy. Taken from [30] . 



total system size. Panel (c) of fig. fT3l shows that the prob- 
abihty of the first group decreases rapidly in the exci- 
tation energy range 2-6 MeV/nucleon, while that of the 
second one increases. The temperatures T associated to 
each class are different, as appears on the related caloric 
curves: the residue-evaporation class shows a Fermi-gas 
behaviour (proportional to T-squared) while the multi- 
fragmentation class is associated to a classical gas (linear 
in T) . The combination of these two behaviours gives rise 
to the plateau zone in the total caloric curve and explains 
the inflexion point of this curve. One is thus dealing with 
a direct bimodal behaviour, with two excitation energies 
associated with one temperature in the transition region. 
This behaviour is an intrinsic feature of the phase space 
population in the SMM. 

3.2 CMD: Classical Molecular Dynamics 

Signals of phase transition were searched for in dynamical 
models. A simple example is a classical molecular dynam- 
ics model with a Lennard- Jones potential implemented 
by Cussol [33]. With such a potential, analogous to the 
van der Waals interaction for fluids, the model includes 
a liquid-gas phase transition. Symmetric collisions of LJ 
droplets with sizes of 50 -|- 50 and 100 -I- 100 are analyzed. 
Systems were prepared in three different conditions: 

" central collisions (small impact parameters), 
— peripheral collisions (all impact parameters but look- 
ing at the forward zone, "quasi-projectiles"), 



— " thermalized" systems (particles are placed in a box of 
volume V/Vq = 8 and released after a time sufficient 
to reach thermal equilibrium). 

Two variables were scrutinized, the size asymmetry be- 
tween the two largest fragments, 77 (eq. |4]), and the mass 
correlation between these same fragments [34]. The re- 
sults, for systems comprising 100 droplets, are presented 
in fig. [14] (top: central collisions, middle: quasi-projectiles 
and bottom: thermalized system). Excitation energies are 
expressed in ESU, ratio between the excitation energy per 
particle and the binding energy of the least bound particle. 

Bimodality - the occurrence of two fragmentation pat- 
terns in a given energy zone - is present in all situations, 
but at different excitation energies: 1 ESU for central col- 
lisions, 1.5 ESU for quasi-projectiles, and ~ 1.8 ESU for 
the thermalized system. It is however worth to mention 
that if the thermalized systems is prepared at higher den- 
sities {p/po = 1 — 1.5) the transition between the fragmen- 
tation patterns also occurs but at lower excitation energy, 
namely < 1 ESU [M]. Cussol attributes the differences 
in the transition energy to the lack of complete thermal- 
ization of any source produced in nuclear collisions, what- 
ever the impact parameter. One can conversely argue that 
this study proves that bimodality is a robust signature of 
phase transition, as it survives even if the system is not 
fully thermalized, although the apparent transition energy 
is displaced. This point will be developed later. 
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ability of events selected on Amax (c) and average temperature 
for these events (d) versus E* for Kr, La, Sn and Au nuclei 
{SMM simulations). Taken from [31]. 
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Fig. 14. Mass correlation between the two largest fragments 
(left) and mass asymmetry, 77, as a function of the excitation 
energy (right), for collisions of LJ droplets. The top panel is 
associated to central collisions; the middle to peripheral ones 
and the bottom panel to " thermalized" systems (see text) [34] . 



3.3 HIPSE: Heavy Ion Phase Space Exploration 

The Heavy Ion Phase Space Exploration model HIPSE 
comprises a full (classical) treatment of the entrance chan- 
nel (nucleus- nucleus potential, AW collisions). It is fol- 
lowed by a random sampling of nucleons in the partic- 
ipant zone from Thomas-Fermi distributions of the two 
colliding nuclei to form fragments in the dense zone |36| . 
Excitation energy is shared among all products, taking 
into account the total energy constraint. Finally a statis- 
tical de-excitation {SIMON code [3S]) of the fragments, 
including QP and QT - if they are still present - is per- 
formed. 

Simulations were done for all impact parameters, to 
mimick a real 50A MeV Xe-|-Sn experiment, then the same 
analysis as in [30| was performed by Lopez et al [37| ; a bi- 
modal structure was observed in the correlation between 
Zmax and the charge asymmetry of the two largest frag- 
ments (eq. 2]) . In a model however one can go further and 
track the origin of the bimodal behaviour: is it due to the 
entrance channel (dynamical effect) or to the de-cxcitation 



step? The first hypothesis was ruled out, as no discontinu- 
ity was found in the evolution of the size of the hot largest 
fragment with the impact parameter: the bimodality was 
clearly attributed to the statistical de-cxcitation of the 
QP. A deeper analysis of the de-excitation stage was then 
achieved through the simulated statistical de-excitation of 
Xenon nuclei of different excitation energies and spins with 
the SIMON code [35]. This is depicted in fig[15l where the 
distributions of the asymmetry variable (eq. [5]) are plot- 
ted for several initial conditions. Increasing the excitation 
energy does decrease the average charge asymmetry, but 
never down to the small values observed in the data. Con- 
versely, if more spin is given to the nucleus, the asymmetry 
variable displays a sharp transition around 60 — 70?i, which 
corresponds indeed to the angular momentum for which 
the symmetric fission barrier vanishes. 

The authors of [37] conclude that, in the HIPSE model, 
the observed bimodality found its origin in the spin rather 
than in the excitation energy transferred to the QP, being 
still a phase transition but not of the liquid-gas type. It 
is worth mentioning that using the SMM model for the 
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Fig. 15. Charge asymmetry distributions resulting from the 
de-excitation of hot Sn nuclei with different initial excitation 
energies (columns) and spin (rows), with the SIMON code. 
Taken from [37]. 



de-excitation stage, the authors also observe bimodality in 
the size of the largest fragment. This is not surprising in 
view of the abovementioned study with the SMM. How- 
ever, this raises the important issue - still under debate - 
of the order parameters (and then the type of the phase 
transition) which govern the bimodality. This point will 
be discussed in the perspectives. 



4 Pending questions 

As seen in the previous sections, bimodality is a very com- 
mon feature in nuclear collisions at intermediate energies. 
It is present in central as well as in peripheral collisions. 
It takes place for a large range of masses, A=40-200. It 
was however mentioned in the course of the text that it 
is experimentally difficult to isolate a source, because of 
dynamical effects leading to a mixture of pre-equilibrium 
products and of QP / QT de-excitation particles. Even if a 
source can be properly defined, one has to verify its degree 
of thcrmalisation. Indeed radial flow was found, particu- 
larly in central collisions, and transparency effects were 
also evidenced [35] ■ It seems however from both experi- 
mental [3D] and theoretical [37] studies that bimodality is 



not mainly driven by dynamical effects. Ambiguities re- 
main in the type of phase transition observed, and conse- 
quently on the definition of a true order parameter. Some 
of these questions were addressed recently and are pre- 
sented in the following. 



Table 1. World-wide experimental results on bimodality in 
July 2005. . 



Results 
from 


Reaction 
centrality 


Source 
size 


Bimodal 
variable 


INDRA 


Central 


~ 200 


Zi - Z2 (eq. [1] H 


INDRA 


Central 


~ 200 


Ziiq - Zgas (eq. [H 


INDRA 


Central 


~ 100 




INDRA 


Peripheral 


160-180 


Zi - Z2 (eq. ID) 


MULTICS/ 
MINIBALL 


Peripheral 


~180 


Zi - Z2 (eq. ID) 


ALADIN 


Peripheral 


~130 


Z-y-Z-2- Zz (eq. O 


NIMROD 


Peripheral 


24-40 


Zliq Zgas 



Table 1 gathers all experimental results on bimodal- 
ity found so far. A glance at the table indicates that bi- 
modality was essentially found in charge asymmetry vari- 
ables comprising the two or three largest fragments of 
each event. Such variables can in some sense be related 
to density difference between a dense (liquid) and a di- 
lute (gas) phase; in some models for instance the Fisher 
droplet model, the largest fragment is assimilated to the 
liquid while all the other form the gas. 



4.1 Are Zmax< ^max< Or the asymmetry order 
parameters? 

Simulations were performed in different frameworks to test 
whether the observables Z^ax, ^max, or the asymmetry, 
reliably sign a phase transition. Let us recall that a bi- 
modality of an order parameter signs the occurrence of a 
first order phase transition in a finite system. Fig. ll6l shows 
the outcomes of three simulations in the transition region 
- when it exists; there is no phase transition in the Ran- 
dom Partitions calculation, while percolation has a second 
order transition and Lattice Gas a first order one |39j . 
The distributions of the largest fragment A^ax evidence 
that Aj-nax only presents a bimodal distribution for the 
canonical Lattice-Gas calculation. This means that A^ax 
is indeed an order parameter of the first order phase tran- 
sition of the Lattice Gas. The distribution presents a wide 
plateau, as expected, in the case of a continuous transition 
(percolation). By contrast, the mass asymmetry, Aasy^ de- 
fined in a similar way as the charge asymmetry (eq.[3|), also 
displays a bimodality (although with a less marked mini- 
mum) for simulations which have a 2"^^ order phase transi- 
tion (percolation, middle column). The conclusion of this 
study is that both A^ax and Aasy clearly signal a phase 
transition - note that none of them presents bimodality in 
a model without phase transition - but Amax is the only 
unambiguous signature of the order of the transition. 
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Fig. 16. Largest fragment A^ax (top) and mass asymme- 
try Aaay (bottom) distributions for three simulations. Left: 
random partitions (no phase transition), middle: percolation 
(2'"*order phase transition) and right: canonical Lattice-Gas 
(1"* order phase transition). Taken from [39j . 



4.2 Order parameters of the liquid-gas phase transition 

If nuclei undergo a liquid-gas type phase transition, then 
the order parameters are known: the energy, the volume. 
In some of the experimental studies cited above, the au- 
thors try to push the analysis beyond the single obser- 
vation of bimodality on the asymmetry variable. As a 
first attempt, in central collisions between Ni and Au at 
52 A MeV [TU], the excitation energies (experimentally 
deduced from the energy balance equation) associated to 
the two fragmentation patterns were found slightly differ- 
ent (by lAMeV) [3D]. This bimodality of the excitation 
energy is an indication in favour of the liquid-gas type of 
the phase transition observed. 

Studies of Au quasi-projectiles were deepened by the 
authors of ref . [3D] : a test of the reliability of the canonical 
picture was accomplished by estimating the apparent tem- 
peratures of the two types of events, from the slope of the 
emitted proton spectra for residue-like events, and from 
double isostope ratios in the multifragmentation regime. 
As seen in fig. [T7] both temperatures are close enough 
in the region where bimodality is present (Etrans = 0.8- 
1.2 A MeV), while the excitation energies, calculated with 
the energy balance equation, arc different. This is expected 
if bimodality has a thermal origin and validates the sorting 
as close to a canonical one. 



4.3 Does bimodality survive out-of-equilibrium effects? 

The influence of non equilibrium effects on signals of phase 
transition was studied in [IT] in the case of incompletely 
relaxed incoming momentum (transparency) and of self- 
similar radial flow. Both effects were indeed recognized in 
experimental data. Fig. [TH] displays results of (canonical) 
Lattice-Gas simulations with different radial flow energies; 

Ap^ /p^ is equivalent to the ratio Sjiowl^th and is var- 
ied from to 1. At the transition temperature, T=0.68e, a 
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Fig. 17. Apparent temperatures of Au quasi-projectiles as 
a function of the normalized transverse energy for residue 
(filled symbols) and multifragmentation events (open symbols). 
Taken from M\. 



bimodality of A„iax is clearly seen in the absence of flow, 
and is still visible even when the flow energy is as impor- 
tant as the thermal energy (top right panel) . The authors 
state thus that radial flow docs disturb the signal, partially 
filling the gap between the two components, but does not 
destroy it as long as the flow does not dominate the global 
energetics. Similar conclusions were drawn in this paper 
in presence of longitudinal flow (transparency effects). 

These two examples illustrate the robustness of bi- 
modality versus external (and realistic) constraints due 
to the dynamics of the collision; similar conclusions can 
also be derived from CMD simulations (see above). 

In experimental data on Au quasi-projectiles [30j . re- 
fined treatments aiming at better isolating quasi-projectiles 
from the mid-rapidity contribution and keeping only events 
where this contribution was smaller were tempted. In all 
cases the bimodal picture comes out better, although it 
occurs for a lower value of the sorting variable (smaller 
dissipation), for a given incident energy. This is again 
an evidence of the robustness of bimodality against non- 
equilibrium effects. 
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Fig. 18. Canonical Lattice-Gas simulations for different tem- 
peratures T around the critical one Tc=0.68e for the distri- 
butions of the largest fragment. Simulations are performed by 
adding an extra radial flow energy Z\p^/p^ between and 1. 
Taken from [2]. 



5 Perspectives 

Bimodality is a very promising signature of first order 
phase transition because of its simplicity and robustness 
against dynamical constraints. It was shown in this con- 
tribution that the signal is quite common for the decay of 
hot nuclei and can be observed in rather different experi- 
mental conditions (central / peripheral collisions, small / 
large source sizes). 

Nevertheless, some open questions need to be answered 
in order to firmly assess the validity of this signal. Several 
strategies can be envisaged in order to progress in this 
direction: 

— cross the observation of the bimodality signal with that 
of all the other proposed signals for the phase transi- 
tion such as critical exponents, scalings (Delta-scaling, 
Fisher scaling, Zipf law), negative heat capacities, or 
space-time correlations (emission times and correlation 
functions). Obviously, when possible, all signals should 
be studied on the same sample of events to minimise bi- 
ases due to sorting. Such cross controls were started by 
the INDRA and NIMROD [22] collaborations. One 
must solve however the problem of the non-equivalence 
of statistical ensembles in some cases. 

— test the effect of sorting. Indeed different ways of sort- 
ing were proposed (impact parameter selectors, com- 
pact shape events, source selection). The robustness of 
any signal will be established if its observation is not 
drastically dependent on the chosen sorting for a given 
centrality for instance. 

— compare the results of different entrance channels for 
nuclear collisions; by using asymmetrical reactions such 



as light ions impinging on heavy targets, or nuclcon/pion- 
nucleus reactions, one may hope to disentangle the dif- 
ferent effects which could possibly govern bimodality. 
By using these very different entrance channels reac- 
tions, the pre-equilibrium / neck contributions can be 
evaluated and even subtracted. Moreover, the effects 
of large collective motions such as radial flow (for cen- 
tral collisions) or spin (angular momentum transfer in 
semi-peripheral reactions) can also be measured. It will 
possibly help to answer to the fundamental question 
of the type of phase transition which is experienced by 
hot nuclei. 



We thank all the nuclear physicists around the world who send 
us their results - published or not. 
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